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1. Introduction 

The spirit that has traditionally guided the construction of approximate density 
functionals has been to collect as much information as possible about the fluid and 
cook up a functional out of it under reasonable simple assumptions. This is the idea 
behind widely used approximations such as the weighted density or the effective liquid 
approximations, which produce a functional out of the Helmholtz free energy and direct 
correlation function of the fluid (see Evans 1992 for a review), or even more recent 
approaches (Zhou 2001a, 2001b) which in addition incorporate the pair correlation 
function. As an alternative to this there stands Rosenfeld's fundamental measure (FM) 
theory, whose aim is to construct the density functional upon geometrical grounds, 
and to obtain as output all the ingredients that classical approaches need as input. 
Since Rosenfeld's pioneering work for hard spheres (Rosenfeld 1989), later extended 
approximately to nonspherical objects (Rosenfeld 1994), this theory has developed far 
from its roots, giving rise to functionals for hard spheres that reproduce extraordinarily 
well the structure of the crystal phase (Tarazona 2000), functionals for the fluid mixture 
of parallel hard cubes (Cuesta and Martinez-Raton 1997a, 1997b) which have been 
applied to study entropic demixing (Martinez-Raton and Cuesta 1998, 1999) and even 
functionals for soft interactions (Schmidt 1999) as well as for nonadditive mixtures 
(Schmidt 2000a, 2000b). 

Rosenfeld's original approach (Rosenfeld 1989), applied to hard spheres, consisted 
of three steps: (i) decomposing the Mayer function into a sum of convolutions of 
one-particle measures; (ii) using these measures to define a set of weighted densities 
and assuming that the density functional depends on the density profile only through 
them, and (iii) determining the functional form by resorting to scaled-particle theory. 
The functional so obtained had an interesting functional structure and worked well 
when applied to problems such as the adsorption profiles at hard walls. However 
it dramatically failed the test of freezing because of some builtin singularities, which 
could be repaired by imposing exact dimensional reduction of the functional to zero- 
dimensional cavities, i.e. cavities which do not hold more than one particle (Rosenfeld 
et al 1996, 1997). This idea was later exploited by Tarazona and Rosenfeld (1997) 
and Tarazona (2000) as a more fundamental idea to construct a functional free of 
divergences. The method amounts to extend to higher dimensions the exact form of 
the one-dimensional hard-rod functional (Percus 1976) and then add the appropriate 
extra terms (with a form dictated by the original FM functional structure) so as to cancel 
whatever exceeds the exact result when the functional is reduced to properly tailored 
zero-dimensional cavities. This method avoids the use of scaled-particle theory while 
at the same time imposes a first-principles constraint (the zero-dimensional reduction) 
to the functional. It is then not surprising that the resulting functional performs much 
better than the original one, for instance in describing the crystal phase (Tarazona 
2000). 

The FM functional for the fluid of parallel hard cubes was first obtained 
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using Rosenfeld's original method (Cuesta 1996) supplemented with the exact zero- 
dimensional reduction requirement (Cuesta and Martinez-Raton 1997a, 1997b). But 
this system seems to be more suitable for FM theory than the hard-sphere one, by 
at least two reasons: first, because the functional turns out to be derivable through 
an extraordinarily simple recipe, and second, because, as we will show for its lattice 
counterpart, the method of cavities leads precisely to that functional. And this remains 
true not only for the monocomponent system but also for a mixture, where the FM 
functional for hard spheres shows its more prominent weaknesses (Cuesta et al 2002). 

In principle nothing precludes us from extending this theory to lattice systems. 
However, when one naively tries to do this immediately runs into difficulties. To 
illustrate them let us suppose we want to follow Ronsenfeld's method. To begin with, 
there is no equivalent to scaled-particle theory for lattice fluids (Soto-Campos et al 
1999) and it is not clear how to devise it, so this makes step (iii) of the derivation 
infeasible. But there is even more: the decomposition of the Mayer function is ambiguous 
in the sense that there is not a unique way to define the fundamental measures. Take, 
for instance, the one-dimensional fluid of hard rods in the continuum. The Fourier 
transform of its Mayer function is 

i i \ sin(oj + aj)q sina^g sina^g 

—Jijiq) = 2 — = 2 cosa^g — h 2 cos a jq , 

q q Q 

where the indices i,j are introduced to distinguish particles and 2a^ is the rod length of 

particle i. As it is clear from this expression, the decomposition suggests itself. For the 

lattice counterpart of this model we have, however, 



-hM) 



sin(aj + aj - \)q 



sin(g/2) 

without being too imaginative, there are several ways of decomposing fij(q), all equally 
plausible: 

t ( x sm («i - |)9 . sin^g x 

-Jij{q) = — . 1 /0 1 cosa^g + . 7— cos(a; - ±)g 



sin(aj — i)g sin(a ; - — ^)g 

= ; — -rA, — cos a,g H ; — —{ — cos a^q 

sin(g/2) 3H sin(g/2) lH 

+ cos a^q cos ajq + sin a^q sin a^g 

sm(ai -\)q sm( aj - \)q 

sm(g/2) 4 sm(g/2) 4 

and many other possibilities. Contrary to the continuum case, this gives no clue as to 

how to define the weighted densities. 

An alternative to this approach is to use the method of cavities, but there is also a 

problem here. The exact one-dimensional functional (Robledo 1980, Robledo and Varea 

1982, Buschle et al 2000a) is very different in its structure from the continuum one, so 

there is no equivalent to defining convolutions of delta-shells with appropriate kernels 

(we will see later that this has no counterpart in a lattice system) . So here again we do 

not have a guide to follow. 
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The method we have employed to achieve the extension and to provide the right 
way to define what could be referred to as lattice fundamental measure theory (LFMT), 
rests on two pillars: first, an appropriate analysis of the extension to mixtures of the 
exact one-dimensional functional of hard rods (which, to our knowledge, we present for 
the first time in this article), and second, a similarity, in the case of cubic lattices, with 
the fluid of parallel hard cubes. After devising the right way to define LFMT we will 
show how to make sense of the theory of cavities for lattice fluids. A further byproduct 
of this theory is that it would allow for an alternative way to define a scaled-particle 
theory for lattice fluids were we able to fill in step (iii) of Rosenfeld's method in such a 
way as to get the same functional. We are currently working along this line. 

The structure of the article goes as follows. In section § we analyze the mixture 
of hard rods on a lattice. This means considering two cases: (a) all rod diameters 
having the same parity, and (b) some having odd and some even diameter. The former 
corresponds to an additive mixture, while the latter is nonadditive. The functional for 
case (a) is derived first, following Vanderlick's et al (1989) derivation for the continuum; 
this method cannot be applied to case (b), which is then obtained through an appropriate 
mapping to an additive mixture in a particular external field. The way in which those 
exact functionals are finally written suggests, by analogy with the fluid of parallel hard 
cubes, the extension to higher dimensions (only in cubic lattices). This is carried out in 
section ^, where it is also proven that this extension has a correct dimensional reduction 
to any lower dimension. It is then show how to extend the method of cavities for these 
lattice models. Finally, in section ^, we present some results from applying the obtained 
functional to some particular two and three-dimensional systems. 

2. Hard— rod lattice fluid revisited 

The exact density functional for the continuum fluid of hard rods was first obtained by 
Percus (1976) and extended to multicomponent mixtures by Vanderlick et al (1989). 
Robledo (1980) and Robledo and Varea (1981) derived the exact functional for the 
lattice fluid of hard rods by means of potential-distribution theory (Widom 1978), but 
they focused on the continuum limit. Buschle et al (2000b) have recently dealt with the 
general problem of the derivation of exact density functionals for one-dimensional lattice 
gases with finite-range pairwise interaction. They made use of a generalized Markov 
property satisfied by the conditional particle distribution probabilities. In this way, they 
rederived (Buschle et al 2000a) the functional for the monocomponent hard-rod lattice 
fluid. 

In this section, we are going to derive the exact density functional of a 
multicomponent lattice fluid of hard rods. The starting point will be the adaptation to 
lattice models of Vanderlick's et al (1989) derivation for the continuum. 

Let us consider a system of hard rods with different sizes, such that the interaction 
between a-type and a'-type rods at positions s and s', respectively, on a one- 
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dimensional lattice (i.e. s, s' G Z) is given by 

{0 if Is — s'\ > cr aa /, 
. ' (2-1) 
oo it \s — s I < cr aa > } 

where a aa ' = \{a a + a a >), cr a G N being the diameter of an a-type rod. Because of the 
particle ordering in one-dimensional systems, the pairwise interaction will appear in the 
grand canonical partition function H via the modified Boltzmann factor, defined by 

e aa ,{s, s') = e- /3 ^' (s ' s ' ) 6(s - s') = G(s - s' - a aa >), (2.2) 

where 0(s) = 1 if s > and otherwise, and (3 = 1/kT with T the temperature and k 
the Boltzmann constant. Note that this eliminates the combinatorial factor 1/N\ of the 
partition function. Let us also suppose that every species a is in chemical equilibrium 
with a particle reservoir of fixed chemical potential p a and that over an a-type particle 
at position s e Z is acting the external potential u a (s). In the same way, the external 
and chemical potentials will only appear in H via 

w aa -{s, s') = e PW - Ua{s)] b~ aa ,5 ss , = w a (s)5 aa ,5 ss ,, (2.3) 

where 8 aa i is the Kronecker delta. With the above definitions we can write the grand 
canonical partition function as 

N-l 



1 + ^Tr 



N=l 



ai( s l) \\ e ociOCi+i( s i, s i+l) w ai +1 ( s i+l) 



"' 

i=l 



(2.4) 



where Tr = ^{(q^ s^} sums over an indices (a, s) appearing in the expression, a« 
denotes the species of the ith particle and Sj G Z its position. Considering e and 
w as operators defined over the (a, s)-space, the term under the sum over iV is just 
Tr { [ u '( eu ') Ar ~ 1 ] Q , a / ( s ! s ')}i thus we can perform the sum and rewrite ( |2.4p as 

S = l + Tr {[w(I-ewr l ] aaf (s,s')}, (2.5) 

/ denoting the identity operator in the (a, s)-space. 

The density profile of species a, p a (s), can be expressed as the functional derivative 
(Hansen and McDonald 1990) 

w a (s) dE 

p a {s) = -^—t, — pr, (2.6) 



which in this case, given the discrete nature of s, is a partial derivative. Using fl2.5|) for 



H in ( p.6|) , we can write p Q (s) as 



A.W = "'"Y"-'" , (2.7) 
where we have introduced the truncated partition functions defined as 

E( J - ew )-'( s ' s '), (2-8) 

(a',s') 

(«',*') 
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(a',s') 



From the above definitions and multiplying by (/ — ew) to the left of S~ (s) and by 
(I — we) to the right of (s), it is straightforward to obtain the following recursive 
relations satisfied by the truncated partition functions 

(2.10) 

S a (S) = 1 + X) ^ ( We )°'«( S '' S )' 

(a',s') 

At this point, our original problem may be posed as to solve (|2.7[) and ( |2.10| ) in such 
a way that w a (s) could be expressed as a functional only depending upon the densities 

p a (s). 

As it was noted by Vanderlick et al (1989) and Percus (1997), the solvability of 
and (|2.10| ) depends on the rank of the Boltzmann factor Q2.2Q regarded as an 



operator in the (a, s)-space. For the additive continuum system it can be shown that 
this is of rank one, and then it is possible to solve those equations. On the contrary, 
for the lattice system we are dealing with, this is not always true. In this case, due 
to the discrete nature of the system, we should distinguish between a mixture of rods 
all of them with even(odd) diameters (additive mixture) and one of rods some of which 
have even diameter and some odd (nonadditive mixture). Indeed, for the former case, 
a aa i is always an integer and, consequently, from (|2.1|) it follows that the corresponding 
mixture is additive. The solution in this case is just a rephrasing of the resolution of 
( p.7| ) and ( |2.10| ) reported for the continuum system (we will sketch it in section |2.1| ). 



For the mixed even-odd case, however, a aa i may take half-integer values. If 
half-integer then \s — s'\ > a aa i in ( [2.1Q is also \s — s'\ > a aa i + 1/2, so effectively the 
interaction is nonadditive. In spite of this, it should be noticed that this nonadditivity 
is peculiar for the lattice system and disappears in the continuum limit. For nonadditive 
mixtures the Boltzmann factor (p.2[) is not a rank-one operator; nevertheless, we will 
show (in section p.2|) how to obtain the exact density functional for that system as a 



particular case of the additive mixture. 



2.1. Additive mixture 

Let us consider a mixture of hard rods such that o a = 2a a + e with o„6l and e = 0, 1 
depending on whether all the rods have even or odd diameters, respectively. In this case, 
the Boltzmann factor ( |2.2| ) is a rank-one operator. If we define the vectors e+(s) = 5 s , aa 
and e~(s) = Q(s — a a — e) and the convolution * as / * g(s) = J2 r f( s ~ r )fi , ( r )> (|2-2j) 
can be written as 

e aa '{s, s') = e~ * e+,(s - s') = e+ * e~,(s - s'). (2.11) 
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e„ is - r) 



{a',s') 



E 

(a',s') 



S+( S ')^( S ')e+( S '-r) 



(2.12) 



We can now use (|2.7|) to eliminate w a '(s') from Q2.12|) ; besides, as e a (s + a a ) = B(s — e), 
hence independent of a, it is handy to define 

^ = ^(^4 (2.13) 
p ± ( S ) = ^p a ( S ±a a ), (2.14) 

a 

which allow us to rewrite equations Q2.12|) in the more suitable form 



p 11 



r=— co 
oo 



(2.15) 



a- r 



From the above expressions it is straightforward to obtain the following system of finite 
difference equations 

oP -(s + l-e) 



AS" s) 



AE + (s) 



S+(s + l-e) 
"H-( S + c)' 



(2.16) 



(2.17) 



with boundary conditions 

S^too) = S, 
-±(±00) = 1, 

where we have introduced the difference operator A, defined by A/(s) = /(s + 1) — /(s). 
This system is exactly solvable and is equivalent to the original recursion relations 
( 2.1U| ). In order to solve fl2.16| ) we have to combine the two equations of the system 
taking into account the discrete Leibnitz rule A(fg)(s) = f(s+ l)Ag(s) + g(s)Af(s) = 
g(s + l)Af(s) + f(s)Ag(s). This yields 

A (S"S + ) (s) — E [p~(s + 1 — e) — p + (s + e)] , 

whose solution, taking into account the boundary conditions Q2.17j ), is 

E-(s)E + (s) = e\i+ J2 [p-(r + l-e)-p + (r + e) 



(2.18) 



(2.19) 



or using the definitions ( |2.14j ), 



s+a a -l+£ 



i - E E pM 

a r=s—a a -\-l—e 



(2.20) 
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This equation permits us to decouple the system ( |2.16| ) and solve it for S ± (s). After a 
little bit of algebra and index manipulation, the solution has the form 



Z + (s) 



11 i_ n (i)( r )- 



(2.21) 



where we have introduced the weighted densities, n^ k \s) (k = 0, 1), defined as 

1 if — a a — k — e < s < a a , 



n^(s) 



otherwise. 



(2.22) 



With (|2.21|) we have basically solved our problem, since from ( 2.21]) and ( 2.17 ) we 
can write the grand potential, f3VL = — InS, as 

'l-nW(s) 



n(°)(s). 

Finally, the free energy density functional, f3J- ', defined as (Evans 1979) 

PF[Pa] = Ptt[Pa] + ^ y^Pa(g) hlW a (s), 



(2.23) 



(2.24) 



takes, from (]27?|), ( CT ) and (ggg) , the form 



X 



s+a a — 1 



s+a Q — 1 



In (1 - n (1) (r)) - In (l - n (0) (r)) . (2.25) 

-s—a a —e r=s— a a — e+1 

For reasons that will be made clear later, we prefer to write this functional in the 
equivalent, but more convenient form, 



(3F\p a \ = (3r d [p a ] + [ $ o {n {1 \s)) - $ (n< >(*))] 



(2.26) 



where {5T lA = X] a 2~^sezPa( s )[l n Pa( s ) — -*-] re f ers to the ideal part of the free energy 
and $o( ? ?) = ?7 + (1 — ^) ln(l — rj) is the zero-dimensional excess free energy for a 
zero-dimensional cavity with mean occupancy rj (Rosenfeld et al 1996). 



2.2. Nonadditive mixture 

As we have mentioned above, the exact density functional for this special system can 
be derived as a particular case of the additive functional ( |2.26| ). Let us suppose that 
we want to study a system of hard rods whose diameters are a a = 2a a + e a , where 
a a e N and e a = 0, 1. Because of the a dependence in e a , the factorization ( |2.11| ) 
cannot be applied to (|2.2j). This difficulty can be overcome if we notice (see figure [I]) 
that this nonadditive mixture can be mapped to an additive mixture of hard rods whose 
diameters are a a = 2o~ a , while their positions are restricted to lay on the even sites of the 
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0=1 
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Figure 1. This figure shows how an additive mixture of hard rods of diameters 2 
and 4 can be made to have the same configurations as a nonadditive mixture of rods 
with diameters 1 and 2. The positions of the upper system are restricted to be on the 
black sites, while the white sites are forbidden. 



lattice (i.e. p a (s) = if s = 2r + 1 with r G Z; this amounts to set an infinite external 
potential on the odd sites, while keeping it arbitrary on the even sites). Therefore, the 
density functional for the nonadditive system will arise by introducing in ( |2.26D the 
density profile 

52p a (r)6 B>2 r, (2.27) 



Pa{S) 



where p a (s) is the density profile corresponding to the original nonadditive system. This 
yields for the weighted densities ( |2.22p , 



n 



(*)/ 



n 



(2s) 



(2.28) 



n[ k) {s)=n (k \2s + l] 

where s G Z, k = 0, 1 and the subindex {e, 0} states that the sum over species in 
the definition ( |2.22 ) is restricted only to those with even or odd diameter, respectively. 
There are thus four instead of two weighted densities for this special case. An alternative 
way of expressing these weighted densities is (s) = ^ a w^ * p a ( s ) = 0, 1), if 
the following nonadditive weights are introduced 

if o a is even, 

(2.29) 

1 + j) if a a is odd. 

With these definitions, the excess free energy density functional for the nonadditive 
mixture (the ideal part remains the same) adopts the form 



w a+*(w)>( a ) 



$ (ni 0) ( s ; 



$0 n, 



(2.30) 



As far as we know, functionals ( 2.26 ) and ( |2.30| ) have never been reported in 
the literature. The additive case is just an exercise of rewriting Vanderlick's et al 
(1989) original derivation for the continuum, but this nonadditive case cannot be 
directly obtained with their method. Besides, even the functional ( 2.26| ) (which for 
the monocomponent case reduces to that obtained several times in the literature) has 
been written in a form that will make of the obtention of an approximate functional in 
higher dimensions an easy task. 
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3. From OD cavities to higher dimensions 



Cuesta and Martinez-Raton (1997a) proved that the one-dimensional excess free energy 
functional in the continuum model can be generated from the zero-dimensional one by 
just applying a differential operator with respect to the particle sizes. The form in which 
the functionals (|2.26|) and (|2.30|) have been written proves that this is also true for the 
discrete case, but now, due to the discrete nature of the system, the differential operator 
becomes a difference operator. This suggests rewriting (|2.26| ) and (|2.30|) as 



E 

b=o,i 



$ 



nf\s] 



(3.1) 



respectively. Here V k is the difference operator V k f(k) = f(l) — /(0). Cuesta and 
Martinez-Raton (1997a) also proved that the FM functional for the <i-dimensional 
system of parallel hard cubes followed the same simple recipe, i.e. it can be generated 
from the zero-dimensional functional by successive applications of a differential operator 
with respect to the particle size along every coordinate axis. The functional so obtained 
satisfies the zero-dimensional reduction requirement. For the equivalent lattice model 
it seems natural to propose a similar recipe. Then, taking into account the continuum 
analogue and ( |3.1| ), our proposal for the FM functional of the <i-dimensional discrete 
parallel hard cube model is 



( 



n 



(k) 



(s)) 



(3.2) 



for the additive case, and 



E 

je{o,i} d 



for the nonadditive one, where k 



(A' 



1: 



nf\s) 



1- 



(3.3) 



,jd) ar e vector indices, 



2\ = Yli=i is the difference operator and the weighted densities are defined as 



n (k) 

(k) 



U. 



(k) 



Pa(S, 



with u>a (s) 



L), 



(fci). 



a 



s^, for the additive system, and as 



for the nonadditive one. 



s ) = Ea w jS * P°( S ) With U £1( S ) 

Functionals ( p.2|) and (|3.3j ) are the main results of this work since they constitute, 
to the best of our knowledge, the first FM functionals for a lattice model. In this 
sense, it can be checked that (i) they depend on the density profile through a set of 
weighted densities defined from one-particle measures which are compatible with the 
decomposition of the Mayer function into a sum of their convolutions; (ii) the virial 
expansion of the direct correlation function is exact up to first order in the density, 
and (iii) they consistently reduce to any lower dimension down to d — 0. Also, it 
must be noticed that their continuum counterpart is recovered in the continuum limit 
of vanishing lattice spacing 5 — > and infinite size a a — > oo with a a 5 = const, (notice 
that the additive-nonadditive differentiation becomes immaterial in this limit). Points 
(i) and (ii) are easy to check, so we will focus in what follows in the crucial point (iii) 
and its connection with an alternative way of obtaining functionals (|3J 



and (|3.3|) 
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a • • o a • • 
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l : : : l I I : 
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Figure 2. Maximal cavities for binary mixtures in d — 1. Large particles can be 
placed at any point, black or white, but small ones can only be placed at black sites. 
The figure shows the additive even (a) and odd (b) cases, as well as the two equivalent 
cavities for one nonadditive case (c) and (d). 

As it was mentioned in the introduction, there is an alternative method to 
Rosenfeld's original one to construct a FM functional. This method will be referred to 
as the zero-dimensional cavities method (Rosenfeld et al 1996, Tarazona and Rosenfeld 
1997, Tarazona 2000) because it is based on the exact dimensional reduction of the 
functional to d — 0. Pursuing the same idea but having in mind the form of functionals 
( p. 261 ) an d (|2.30|) , we have devised a similar method suitable for lattice models. The 



method amounts to (a) directly extending the exact zero-dimensional functional to 
higher dimensions (note here a difference with the continuum case (Tarazona and 
Rosenfeld 1997), where it is the one-dimensional functional what is extended to higher 
dimensions), and then (b) adding the appropriate extra terms (in the way suggested by 
( [2.261 ) an d Q2.30|) ) in order for the functional to recover the exact form when applied to 
zero-dimensional cavities. 

To illustrate this procedure we must first recall that we are dealing with a mixture, 
sometimes nonadditive, and, contrary to the monocomponent case, it is not very clear 
what a zero-dimensional cavity means for a mixture. So let us start by defining what 
we will understand by that. In the monocomponent case, a zero-dimensional cavity 
is simply a set of connected points in the (^-dimensional lattice such that if a particle 
(its center of mass) is placed at one of them no other particle can be placed at any 
other point of this set. For one-dimensional hard rods of length 2a, the largest zero- 
dimensional cavity ('maximal' cavity) is the interval [—a + l,a], whereas for hard rods 
of length 2a + 1 it is [—a + l,a + 1]. In d dimensions the maximal cavity is the set 
[—a + 1, a] d or [—a + 1, a + l] d , depending on the length parity. 

For a mixture, a zero-dimensional cavity can be defined as a collection of sets, one 
for each species, such that if a particle of any species occupies one of the points of its 
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corresponding set, no other particle of the same or different species can be placed at any 
point of its corresponding set (figure |2] will help to clarify this concept). In d = 1 and in 
the additive case, this means the collection of sets {[— a a +l, a a ]} a , if 2a a is the rod length 
of species a (figure 0a), or {[— a a + 1, a a + l]} a , if this rod length is 2a a + 1 (figure 0b); in 
the nonadditive case, however, indexing by a the species of even length and by a' those 
of odd length, there are two maximal cavities: {[— a a + 1, a a ], [— a a /, a a r]} a>a i (figure 0c) 
and {[— a a + 1, a a ], [—a a / + 1, a a i + l)} a ,a' (figure ||d). In d dimensions the cavities for 
the additive case generalize as in the monocomponent case, while in the nonadditive 
case there will appear 2 d different cavities, depending on which of the two choices above 
we select for each dimension. Because of this proliferation of maximal cavities in the 
nonadditive case, it is convenient to discuss both types of mixtures separately. 



3.1. Additive mixture 

For density profiles constrained to be zero outside the points of a maximal zero- 
dimensional cavity it is easy to check (by applying the definition) that the weighted 
densities ( |2.22| ) satisfy the relationships 

n«(«)=» (0) (« + l) ( S <0), 
n «( s ) = n (0) (s) (s > 0). 

These are the keystone of the method. As explained, the method starts off from the 
straightforward extension of the zero-dimensional excess functional to d = 1, ^q 1 ^. As 
n^(0) = r), the occupancy probability of the cavity, 

s o M M = £M» (1) 00)- (3.5) 

sGZ 

Now we apply this functional to a zero-dimensional cavity and, by means of ( |3.4| ) rewrite 
it as 

< } [p] = $ (n«(0)) + (n(°)(5 + 1)) + (n^(s)) 

s<0 s>0 

= $o(n«(0))+^$o(n( )( S )). (3.6) 

But $ (n (1) (0)) at the r.h.s. of O is the exact zero-dimensional excess free energy for 
the cavity, so the method dictates that in order to obtain the excess functional for the 
one-dimensional system we must substract to the form we started from all the spurious 
terms. This leads to 

which coincides with (12.261). 



Let us now illustrate how the method works for higher dimensions as well. For the 
sake of simplicity we will restrict ourselves to do it only for d = 2. The corresponding 
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relationships between the weighted densities of density profiles vanishing outside a zero- 
dimensional cavity are 

n^ 1 \s 1 ,s 2 )=n^ \s 1 ,s 2 + l) (s 2 < 0), 

n^(s 1 ,s 2 )=n^°\s 1 ,s 2 ) (s 2 > 0), 

with k — 0,1, and similar ones if the coordinate axes are interchanged. Using the same 
notation as for the one-dimensional case and repeating the steps in (|3.6|) for s 2 we get 

*?*\p] = 5>° (™ (1,1) ( s )) = E $ ° + E $ o (^ (1 ' 0) ( s )) • ( 3 - 9 ) 

sgZ 2 siGZ seZ 2 

Repeating the trick with Xlsez 2 (n^^^s)) we get 

= $ o (n (0,1) (s)) - E $ o (^(^ °)) - E $ o (^ ( °' 0) (s)) • (3.10) 

sgz 2 siez sgz 2 

Adding (|0| to ([p]) and noticing that 

£ $ (n^fo, 0)) - J2 $ o (n^V. 0)) = $ (n^(0)) (3.11) 
yields 

^ M) [p] = $o (nW)(0)) 

+ E [ $ o (* (1,O) 00) + $ o K^)) " $ o (^ ( °' 0) («))] • (3-12) 

sGZ 2 

Again, $o (^^ 1,1 - ) (0)) is the exact zero-dimensional excess functional, therefore 

/^ cx [p] = E[ $ ° (" (1,1) 00) 
sez 2 

-$ (n^)(s)) - $ (n^Cs)) + $ (n^(s))] , (3.13) 

which coincides with ( |3.2j ) for d = 2. This provides the rationale for ( |3.2j ) for arbitrary 
dimension, and proves in passing that functional ( |3.2| ) has a consistent dimensional 
reduction to any lower dimension. 

3.2. Nonadditive mixture 

In the additive case it is always the weighted density n^(0) what gives the occupancy 
probability, n, of the (unique) maximal zero-dimensional cavity, so extending the zero- 
dimensional excess functional is as straightforward as equation ( |3.5| ) shows. In the 
nonadditive case, however, this is not as simple. Due to the degeneracy of maximal 
cavities, which weighted density yields the occupancy probability depends on which 
cavity we choose. So before proceeding with the construction of the functional we need 
a unique expression in terms of the one-dimensional weighted densities that provides 
the excess free energy of any maximal cavity. One such expression is 

$Hp] = *o (^(0)) + $o (4 1} (0)) - $o (nf\0)) , (3.14) 

In order to make it easy to understand why this is so and the discussions to come, 
we have introduced a diagrammatic notation for a particular case in figure We will 
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(a) 



GOQ 
W (D 




OOO GO COO 



w (1) 

vv o 



=n?(0) 
=n«:(0) 
=n?(0) 



w<?) 



w< ' 
VB o 



(c) 



' ' @®(^ - =n ( 1 1) (0) 

' OS® ^ = n o , (0) 
i i ®® ^ =n?(0) 



Figure 3. (a) Diagrammatic representation of the weights corresponding to a binary 
mixture of rods of lengths o\ = 3 and o<i = 2. Circles window points whose density 
contribute to the weighted density: empty circles for the large rods and shaded circles 
for the small rods, (b) and (c) Result of convoluting the weights with the density 
profiles of zero-dimensional maximal cavities. 



consider a binary mixture of rods of lengths a% = 3 and a 2 = 2, which is prototypical 
of this case. For this mixture we represent the weights as a chain of empty and shaded 
circles (see figure [3]a). The effect of convoluting a weight with a density profile is to 
select the density at those sites overlapped by the empty circles for the large rods and 
by the shaded circles for the small rods, and add them up. The result of this operation 
is the corresponding weighted density. The application of these weights to each maximal 
cavity is depicted in figure ||b and c, respectively. Figure ||b illustrates that n^(0) = n 
and n^(0) = nf\o), so (|3.14|) reduces to &o(r]), as expected. On the other hand, 
figure ||c illustrates that n^\o) = n and n^\o) = n^O), so again ( [3.14|) reduces to 

We can now carry on with the method as in the additive case and propose 

® { o d) [p\ = £ [*o (n?(s)) + $ (n«00) - $ (n?\s))] . (3.15) 



Whichever the maximal cavity we choose, from definitions ( |2.28| ) and the 



relationships (^T4|) the following analogous relationships between the nonadditive 
weighted densities hold 

4\s)=nf\s), n?(s) = n$\s + l) (s < 0), 

njp(s) = 4 0) (s), n£\s) = nf\s) (s > 0). 
Proceeding as we did for the additive case, we use the above relationships to rewrite 
( CT ) as 

s<0 ' s>0 

= *rW + £*o(4 0) (*)), (3.17) 
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Figure 4. (a) Diagrammatic representation of the weights corresponding to a binary 
mixture of squares of edge lengths ax = 3 and 02 = 2 (see caption of figure [5]a to 
understand the circle notation), (b) Maximal cavities for that system. Note that the 
result of convoluting the weights in (a) with the density profiles of (b), n^ k \o), is 
obtained by superposing the lower leftmost shadow circle of the weight with the lower 
leftmost black point of the cavity. 



where we have made use of ( 3.14Q , which gives the zero-dimensional excess functional; 
so the method gives 



P^[p a ]=J2\^(n?(s))+%(nS\s 



n 



(0), 



$0 (n, 



(3.18) 



which coincides with ( |2.30| ). 

In order to extend the cavities method to dimensions higher than one, we have to 
generalize expression ( |3.14 ). Due to the cubic symmetry of our system this can be done 
just by applying the same rule on each coordinate axis. For the sake of simplicity, we can 
define the operator T£ as T£f(k,j) = /(l, 1) + /(1,0) — /(0, 1) and its d-dimensional 



version as 7^ =Yll 



=1 1 h i ■ 



Then the natural generalization of (|3.14Q is simply 



T^„(»< k) (o; 



and the extension of the zero-dimensional excess functional to d dimensions 



(3.19) 



(3.20) 



Again, as a representative of a higher dimension, we will apply the cavities method for 
d — 2. In figure § it is shown the diagrammatic representation of the weights for a 
binary mixture of squares of edge lengths o\ — 3 and 02 = 2, as well as all maximal 
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cavities for that system. Expression ( |3.19| ) for d = 2 accounts for all cavities in figure |]b. 
For example, if we take the cavity at the top the application of the weights in figure 0|a 



gives ngfi^O) = V , rJff'{0) = n^O), n\^>(0) = ^'(0), n^>(0) = <f(0) and 



,(i.i)/ 



,(o,o) 



,(0,1). 



,(1.1), 



,(1.°)/ 



n i i (0) =n ii (0), so (|3.19| ) reduces to <&o(r)). By symmetry the same happens for any 
of the cavities. 

Proceeding with the method, from ( |3.20| ) we will start with 
= E [ $ o + *o (ngfts)) - % 



sez 2 



+$ (nSf(s)) + $ (< 1} (s)) - % (4!i 0) (s; 



(o,i). 



,(0,0), 



Now, the relationships (|3.16|) become 



n 



^(sx, s 2 ) = n^i'isx, s 2 ), n^{ l> (sx, s 2 ) = n^'(sx, s 2 + 1) (s 2 < 0) 



,(fc,o). 



,(fe,i) 



(fc,0). 



(3.21) 



(3.22) 



njfd ( s i> s 2) = n^ 0) (si, s 2 ), n^T^i, s 2 ) = n^\s u s 2 ) (s 2 > 0), 
and those obtained by interchanging the coordinate axes. As before, by making use of 
these relationships in ( |3.21| ) we get 

CV] = £^o (nf\sx,0) 



+ E [*° (4:o 0) (s)) + $ Uf(s)) - $ (n?f (s) 



(3.23) 



sez 2 



where it must be noticed that the second term in the r.h.s. of (|3.23| ) is just 



n% Kf (s) 



By the same argument, we have 



,(o,o). 



(3.24) 



= E T ^ (<f&) ~ E (n$f(sx,0)) - E$o 

sez 2 siez sez 

Adding ( p.24|) to ( |3.23|) and noticing that because of the zero-dimensional reduction of 
the one-dimensional case proved above 



E (»«( 8ll 0)) «f\sx,0))} = $rV] 



(3.25) 



we obtain 



+Eh«i 1) ( s ))+^Ko 1) ( S 



sez 



EKRo 0) (s; 



$o(n£ 0) (s))+$ (^7(s 



,(0,0), 



(3.26) 



Thus, the cavities method produces 



sez 2 
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(3.27) 



which is (3J3) for d = 2. 

The extension of this argument to arbitrary dimension proves that (|3.3|) is the 
functional produced by the cavities method. In addition, it is also proved its consistent 
dimensional reduction to any lower dimension. 



4. Some applications 



In order to show the applicability of the theory exposed in the previous sections, in 
this one we present the results obtained for some particular two and three-dimensional 
systems. 



4-1. Hard square lattice gas 

The hard square lattice gas is just the lattice gas with first and second nearest-neighbour 
exclusion in the simple square lattice. This system has been widely studied (Bellemans 
and Nigam 1967, Ree and Chesnut 1967, Nisbet and Farquhar 1974, Slotte 1983) but no 
definitive conclusion has been reached about its phase behaviour. The results reported 
cover the whole spectrum, depending on which theory has been used to study it. Some 
authors have claimed it to have a second order transition (Bellemans and Nigam 1967, 
Slotte 1982), others a weaker (third order) transition (Bellemans and Nigam 1967, Ree 
and Chesnut 1967) and others even no transition at all (Bellemans and Nigam 1967, 
Nisbet and Farquhar 1974, Slotte 1982). Nowadays, the only thing that seems clear is 
the structure of the ordered phase: periodic along one coordinate axis while uniform 
along the other (i.e. columnar). However, no simulations have ever been performed in 
order to clarify the disputed order of the transition. 

The result of FM theory is shown in figure [|. The system is found to have a second 
order transition from a fluid phase to a columnar phase. The value of the packing 
fraction at the transition is r/ c = 3 — = 0.764 and the chemical potential at that 
point (/3fi) c = 2.41. This result is compatible with that obtained by Bellemans and 
Nigam (1967) using Rushbrooke and Scoins's (1955) method. They also found a second 
order transition, but at r] c = 0.807 and (/?/i) c = 2.85. The equation of state they 
obtained with that method is compared with the FM result in figure [5[ There is good 
agreement between both theories, being only in the critical region where the two curves 
deviate from each other. 

It is worth mentioning that within the FM approximation this system is analytically 
solvable. Omitting the calculations and denoting L(x) = xlnx, we have obtained for 
the Helmholtz free energy density of the fluid branch 

P$(ri) = L(t]/A) + L{\ - »j/4) - 2L(1 - rj/2) + L(l - V ), (4.1) 
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Figure 5. Equation of state (reduced pressure, f3P, vs. packing fraction, 77) from 
FM theory (solid line) and from Rushbrooke and Scoins's method (dashed line). 



and for the columnar branch 

P$(ri) = l - [L( Vl /2) + L(l - Vl /2) + L( V /2 - Vl /2) + L(l - n/2 + m /2) 

-2L(1 - 77/2) -L(l- V + 771) + 2L(1 - 77)] , 
where the occupancy probability of a sublattice is given by 



771(77) 




8n 2 + I677 - 8 



(4.2) 



(4.3) 



2 V 4n 

Clearly more work, either by means of simulations or by other methods, is needed 
in order to resolve the nature of the transition. This notwithstanding, we do not expect 
good agreement except at low or high pressures; critical points are difficult (sometimes 
impossible) to capture, even approximately, with mean-field-like theories like this one. 
Nevertheless, we think of the FM approach to be valuable in the sense that (i) it is 
possible to obtain an analytic approximation; (ii) the result obtained is compatible 
with one previously reported using a different theory, and (iii) it correctly predicts the 
structure of the ordered phase. 



4-2. Monocomponent and bicomponent hard cube lattice gas 

We have also applied the theory to three-dimensional systems. In this section we present 
the results obtained for the monocomponent systems of parallel hard cubes of edge 
lengths a = 2 and a = 6, as well as for the binary mixture of them. These being three- 
dimensional systems, the ordered phases that compete in the phase diagrams are three: 
smectic (one-dimensional order), columnar (two-dimensional order) and solid (three- 
dimensional order). The period of the ordering is determined from the wave-vector at 
which the first divergence of the structure factor occurs. In the first case this is n, hence 




Figure 6. Equation of state obtained from FM theory for the monocomponcnt hard 
cube lattice gas with edge lengths a = 2 (a) and a = 6 (b). Phases are labeled F 
(fluid), Sm Q (smectic) and C a (columnar), where a — 2,6,7 stands for the ordering 
period. The inset in part (a) shows a very narrow first order Sm2-C2 transition. 



the period is 2. As this is also the period of the closest-packed phase it is the only one 
that is found in stable phases. In the second case the wave-vector is 2tt/7, implying a 
period 7; however, in this case the closest-packed phase has period 6, so both, period-6 
and period-7 phases must be considered. The same holds for the binary mixture. 

The results for the monocomponent systems are shown in figures |]a {a = 2) and 
||b (<t = 6). The first system undergoes, at rj c = 0.568, a second order transition from 
a fluid phase (F) to a period-2 smectic one (Sm.2). At a higher value of the packing 
fraction the system is found to have a very narrow first order phase transition, where 
the smectic phase at ?7sm = 0.673 coexists with a period-2 columnar phase (C2) at 
?7c = 0.677 (see the inset in figure |6|a). For the a = 6 system we find a more complex 
scenario with several ordered phases, due to the proximity of their free energies. There 
are six ordered phases in competition: smectic, columnar and solid with both periods 6 
and 7. Only three of them, the two columnar and the period-7 solid, are stable in some 
region. The other three are just metastable. Again, the columnar phase is the most 
stable at high density (notice the difference of the high density ordered phase of these 
lattice systems with respect to that of the continuum model (Groh and Mulder 2001)). 
Upon increasing density, the sequence of transitions is as follows: first, a second order 
one at r] c = 0.402 from a fluid phase (F) to a period-7 solid phase (S7); then a first 
order transition in which the S7 solid coexists, at i]s 7 = 0.617, with a period-7 columnar 
phase (C7), at f]c 7 = 0.631; finally, another first order transition with the C7 columnar 
at T]c T = 0.656 coexisting with a period-6 columnar phase (Cq) at r/c 6 = 0.827. This 
latter phase extends up to the closest-packing. 

The competition of different ordered phases in the latter system is remarkable, but 
it produces a real mess when the two kind of cubes are mixed. The phase diagram for this 
binary mixture (pressure vs. composition) is shown in figure [j]. All possible phases are 
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Figure 7. Phase diagram (reduced pressure, j3P, vs. composition, x) of the lattice 
binary mixture of hard cubes (size ratio 6:2). Composition is defined by cc = T)l/t)i 
where rj — i% + f]s is the total packing fraction of the large (L) and small (S) cubes. 
Phases are labeled as in figure ^. The dotted line corresponds to the spinodal of the 
uniform fluid. For 0.81 < x it marks a stable continuous F-S7 phase transition. 

present there. A few particularly complex regions are shown in more detail. Remarkable 
features of this phase diagram are the very wide F-C6 phase separation dominating the 
lower part of the diagram, and the extraordinarily narrow S6-Cg coexistence, extremely 
close to the monocomponent fluid of large cubes. The latter tells us about the proximity 
of the free energy of the solid phase to that of the columnar in this monocomponent 
system (a small amount of 'impurities' breaks the translational symmetry along the 
columns of the columnar phase transforming it into a solid). The former has a very 
important consequence concerning the demixing scenario for this system. 

This binary mixture has been reported, from simulations (Dijkstra and Frenkel 
1994), to demix into two fluid phases. As a matter of fact, this is, as far as we know, the 
only existing evidence of an entropy-driven fluid-fluid demixing in an additive binary 
mixture. In all other systems of the same kind, e.g. hard spheres (Dijkstra et al 1999) or 
parallel hard cubes in the continuum (Martinez-Raton and Cuesta 1998, 1999), fluid- 
fluid demixing is preempted by freezing of (at least) the largest component. Our results 
contradict this simulations and fit in what seems to be the standard scenario of demixing 
for additive binary mixtures, namely fluid-ordered phase demixing (the fact that the 
ordered phase is a columnar is just a peculiarity of the model). A more detailed analysis 
of this result can be found in Lafuente and Cuesta (2002). 
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5. Conclusions 

In this article we have proposed a formulation of FM theory for lattice models. As it 
was mentioned in the introduction, it is not possible just to carry on with Rosenfeld's 
original formulation but in discrete space. This is so because of the lack of a scaled- 
particle theory for lattice models. In fact, FM theory in the continuum can be regarded 
as a generalization of scaled-particle theory to inhomogeneous phases. For this reason, 
we believe that the knowledge of a lattice FM theory could guide the derivation of a 
lattice scaled-particle theory. We are presently working in that direction. 

In the recent reformulations of FM theory the scaled-particle equation of state is 
no more an input but an output of the theory. Two fundamental ideas appear as the 
keystone of these formulations: (i) the exact density functional must have a consistent 
dimensional reduction to any lower dimension, and (ii) the exact zero-dimensional limit 
is necessary for the correct description of the freezing transition. The latter turns out 
to be a very stringent constraint for the functional and seems to contain most of the 
information needed for the construction of a higher dimensional functional. In this 
context, the zero-dimensional cavities method appears as a scheme for the construction 
of FM functionals. The lattice FM theory we have proposed is the generalization of 
this method suggested by the exact form of the functional of the one-dimensional 
mixture of hard rods on a lattice. Although this system has been previously studied 
in the literature, the originality of our derivation is that (i) it is the first time that the 
distinction between additive and nonadditive mixtures have been noticed and resolved, 
and (ii) the exact functional is written in a form in which the zero-dimensional functional 
arises naturally. The cavities method can be extended to other particle shapes and other 
lattice structures. This will be the content of a future work. 

We have applied the functional to two and three-dimensional models and found 
very reasonable results. In two dimensions these are comparable to those previously 
obtained by other means, and in three dimensions they agree with previous simulations 
data, in the case of the mixture (Lafuente and Cuesta 2002), although with an utterly 
different interpretation which has yet to be confirmed. One remarkable feature of all 
three-dimensional models is the very rich and complex phase diagrams they possess 
(specially the mixture), a complexity which the present theory has nevertheless allowed 
us to tackle without problems. 
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